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A method based on Maximum-Entropy (ME) principle to infer photon distribution from 
\ on/off measurements performed with few and low values of quantum efficiency is addressed. 

. The method consists of two steps: at first some moments of the photon distribution are 



retrieved from on/off statistics using Maximum-Likelihood estimation, then ME principle 
is applied to infer the quantum state and, in turn, the photon distribution. Results from 
00 \ simulated experiments on coherent and number states are presented. 



> ■ I. INTRODUCTION 

o ■ 
o 

Besides fundamental interest, the reconstruction of the photon distribution of an optical signal 



p, plays a major role in high-rate quantum communication schemes 
is required for implementations of linear-optics quantum computing 



msed on light beams l], and 



o 
o 

CL|| have been indeed developed, though their current operating conditions are still extreme At 



Effective photon counters 



■ present, the most convenient method to infer photon distribution is perhaps quantum tomography 



J], which have been applied to several quantum states P, |^ with reliable statistics. However, the 
tomography of a state needs the implementation of homodyne detection, which in turn requires 
^ , the appropriate mode matching of the signal with a suitable local oscillator at a beam splitter. 

■ As a matter of fact, quantum tomography has been developed to gain a complete characterization 

of the signal under investigation, and may not be suited in case where we are interested only in 
obtaining partial information, as for example the photon distribution (|^. 

An alternative approach, based on an array of avalanche on/off photodetectors with different 
quantum efficiencies has been suggested ^ and demonstrated with fiber-loop detectors In this 
scheme, repeated preparations of the signal are revealed at different quantum efficiencies, and the 
resulting on/off statistics is then used to reconstruct the photon distribution through Maximum- 
jikelihood estimation. The statistical reliability of the method has been also analyzed in details 

In this paper we want to further reduce the experimental requirements. We assume that 
avalanche photodetection may be performed only with low values of the quantum efficiency and, in 



2 



addition, that only few of those values are available. Then we analyze whether the photon distri- 
bution may be inferred from this scarce piece of information. We found that the use of Maximum 
Entropy principle, together with Maximum Likelihood estimation of moments of the distribution, 
provides an effective method to obtain the full photon distribution. 

In Section^ we describe in details the two-step MaxLik-MaxEnt method, whereas in Section 
mil results from numerically simulated experiments on coherent and number states are reported. 
Section ITVl closes the paper with some concluding remarks. 

II. PHOTON DISTRIBUTION USING LOW EFFICIENCY AVALANCHE DETECTORS 

Maximum Entropy (ME) principle is a powerful tool to infer the density matrix of a quantum 
state when only a partial knowledge about the system has been obtained from measurements. Let 
us call O the set of independent operators Oj/ summarizing the measurements performed on a 
given systems. O is called the observation level gained from the measurements. The ME 
principle states that the density matrix compatible with a given O is that maximizing the Von 
Neumann entropy while reproducing the known values of the operators in C In order to enforce 
this condition the inferred density operator will depend on a number of parameters whose values 
must be properly determined. ME principle is a way of filling our ignorance about the system, 
without assuming more than what we have learned from measurements. The actual form of the 
inferred state heavily depends on the observation level O: it can range from a (nearly) complete 
knowledge of the state, as it happens for quantum tomography, to the (nearly) complete ignorance. 
An example of the latter case is the knowledge of the mean photon number alone, for which ME 
principle yields a thermal state. In general, the ME density operator qme that estimates the signal 
p for a given observation level is given by 



where Z = Tr [exp {— AjyOi/}] is the partition function, and the coefficients Xi, are Lagrange 
multipliers, to be determined by the constraints 



where (O;^) are the expectation values obtained from the measurements. Suppose now that we 
would like to measure the photon statistics of a given signal p using only avalanche photodetectors 
with efficiencies t]^, v = 1,...,N, i.e. we perform on/off measurements on the signal with N 




(1) 



Tr [qme Oy] 



dx^ logZ = (0,) 
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different values of the quantum efficiency. The statistics of the measurements is described by the 
probability of the off events when the quantum efficiency is rji, 

Pv = Tr[/)n,,] 

oo 

= Y.{1 - tj^Tp^ . (2) 

n=0 

In Eq. (O Pn = ('^IpI'^) and liy = HosiVv), where the probability measure (POVM) of the 
measurements is given by 



nofr(r?.) = X;(l-r/,r |n)(r 



.- , , , {n\ 

71=0 

^oniVu) = I-Uoffiv.) ■ (3) 

From the ME principle we know that the best state we can infer is given by that, in this 
case, reads as follows 

QME = exp < - y X^U^ > . (4) 



N >| 



Explicit equations for the \y are obtained expanding the above formulas in the Fock basis; we have 

Er=o(i - ^^T exp {- (1 - ^m)"} 

Vv = -( ^— s ■ (5) 

Eq. © can be solved numerically in order to determine the coefficients Aj^ and in turn the ME 
density operator. In the following we consider situations where the experimental capabilities are 
limited. We suppose that on/off measurements can be taken only at few and low values of the 
quantum efficiencies. In this case the statistics (j2| can be expanded as 



n 

1 



^n{n-l){n-2)rjlpn + ...} ■ (6) 



Summing the series we have 

1 



py = l-r]yNi + ^{N2- Ni)r]l + 
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-{N3-3N2 + 2Ni)r], , ... 
b 

1 - Niir], + ^r?2 + ^rjl + ...)+ N2{^vl + + • • •) 

-Nsi^rjl + ...) + ... (7) 
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where 

Nk=^ = Tt\p (a^a)''] , (8) 

are moments of the photon distribution. By inversion of Eq. JTI, upon a suitable truncation, we 
retrieve the first moments of the distribution from the on/ off statistics at low quantum efficiency. 

The most effective technique to achieve the inversion of the above formula is Maximum- 
Likelihood. The Likelihood of the on/off measurement is given by 

N 

L = \{p'i-^{l-p.f"-^% (9) 
i/=i 

where riy and Nu are respectively the number of off events, and the total number of measurement 
when using a detector with efficiency r]^. In practical calculations it is more convenient to use the 
logarithm of ©, that reads as follows: 

C = logL 

N 

= ^ {riy logp^ + (A/; - n^) log (1 - p^)) . (10) 
Without loss of generality we can set My = AA, and divide (fTU]) by M , obtaining: 



where f^ are the experimental frequencies of the off events. Substituting Q in we find an 
expression for the renormalized likelihood as a function of the moments N^. 

Maximization of Ulljl over the parameters N}^ leads to the following set of optimization equations: 

1^ = vf/^:_illA^^ = o (12) 

If we stop the expansion in ((T)) at the second order (see below) the derivatives in H12|) take the 
form: 



N 



(13) 



dN2 2 

The system (|12() can be easily solved numerically given the on/ off statistics as well as the number 
and the values of the quantum efficiencies used during the experiments. 
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After having determined the first moments, the density matrix of the signal, according to 
Maximum Entropy principle, is given by 



QME = Z ^exp |-^A,,(a'fa)''| 



(14) 



with Ajy to be determined as to satisfy 



Nk 



Tr 



QME (a'^a)^ 



(15) 



E^=oexp {- ELi 

Notice that the unknowns Xi, are contained also in the denominator and that a suitable truncation 
should be adopted (which can be easily determined by imposing normalization on the ME density 
matrix). 



III. SIMULATED EXPERIMENTS 



We have performed simulated experiments of the whole procedure on coherent and number 
states. For this kind of signals the first two moments are sufficient to obtain a good reconstruction 
via ME principle because their Mandel parameter Q = —1 + (A'^2 — ^D/^i is less than or equal 
to the average photon number Ni [l^ . while squeezed states, for which Q > Ni, are ruled out, 
requiring the knowledge of a considerably larger set of moments (in principle, all of them). A 
number = 5 of values of the quantum efficiency ranging from 1% to 5% have been enough in 
order to achieve good reconstructions. Our results are summarized in Figs. Notice that a 

faithful photon statistics retrieval needs a sufficiently accurate knowledge of the Pu's: our results 
are obtained using M = 10^ observations for each rj^. This condition can be relaxed if we increase 
the number A^ of probabilities measured, but their range of values should be kept narrow because 
Eq. must hold. 

Finally in Fig. |l]we check the robustness of the ME inference, against errors in the knowledge of 
the parameters A^i and N2 that may come from ML estimation in the first step. The quality of the 
reconstruction has been assessed through fidelity 

n 

between the inferred qn = {n-lgMEln) and the true p„ = {n\p\n) photon distribution. As it is 
apparent from the plot the reconstruction's fidelity remains large if the relative errors on both 
parameters are about ±5%. 
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FIG. 1: Reconstruction of the photon distribution of a coherent state with |ap — 1 (black), and comparison 
with the theoretical values (light gray). Simulated on/off measurements have been performed with N = 5 
values of 77 ranging between 1% and 5%. The number of experimental measures is 10^ for each rj. Fidelity 
of the reconstruction is larger than F = 99%. 



P(n) 




FIG. 2: Reconstruction of the photon distribution of a coherent state with |ap = 2 (black) and comparison 
with the theoretical values (light gray). Simulations parameters are the same as in Fig. ^ Fidelity of the 
reconstruction is larger than F — 99%. 



IV. CONCLUSIONS 



In this paper we have shown that the photon distribution of a light signal p, with Mandel 
parameter lower than or equal to its average photon number, can be reconstructed using few 
measurements collected by a low efficiency avalanche photodetector. The on/off statistics is used 
in a two steps algorithm, consisting in retrieving the first two moments of the photon distribution 
via a Maximum-Likelihood estimation, and than inferring the diagonal entries of p using of the 
Maximum Entropy principle. The last step imphes the solution of a nonlinear equation in order to 
shape the statistic to reproduce exactly the moments obtained in the first estimation. Though this 
last process may be delicate, we showed with simulated experiments that it yields sound results 
when applied to coherent and number states. Finally we demonstrated that the method exhibits 
a sufficient robustness against errors deriving from the Maximum-Likelihood estimation. 
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FIG. 3: Reconstruction of the photon distribution of the number state |2) (black) and comparison with 
the theoretical values (light gray). Simulations parameters are the same as in Fig. Fidelity of the 
reconstruction is = 98.5%. 



10.8 11.4 12.0 12.6 13.2 




3.0 



N2 



FIG. 4: Fidelity of the photon distribution reconstruction for a coherent state jap = 3, as a function of the 
average photon number A^i and its second moment iV2. The true values for the two parameters are printed 
in italics. Simulations parameters are the same as in Fig. ^ Notice that the first two entries in the top 
line are set to because, in these cases, the Mandel parameter q is smaller than -1, so they are not physical 
states. 
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